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ABSTRACT 



Context. Were a potentially habitable planet to be discovered, the next step would be the search for an atmosphere and its 
characterization. Eventually, surface conditions, hence habitability, and biomarkers as indicators for life would be assessed. 
Aims. The Super-Earth candidate Gliese (GL) 581 d is the first potentially habitable extrasolar planet discovered so far. 
Qh ' Therefore, GL 581 d is used to illustrate a hypothetical detailed spectroscopic characterization of such planets. 

£jj , Methods. Atmospheric profiles from a wide range of possible ID radiative-convective model scenarios of GL 581 d were used 

to calculate high-resolution synthetic emission and transmission spectra. Atmospheres were assumed to be composed of N 2 , 
CO2 and FLO. From the spectra, signal-to-noise ratios (SNR) were calculated for a telescope such as the planned James Webb 
Space Telescope (JWST). Exposure times were set equal to the duration of one transit. 
q , Results. The presence of the model atmospheres could be clearly inferred from the calculated synthetic spectra due to strong 

5^ ■ water and carbon dioxide absorption bands. Surface temperatures could be inferred for model scenarios with optically thin 

spectral windows. Dense, C02-rich (potentially habitable) scenarios did not allow for the characterization of surface tem- 
peratures and to assess habitability. Degeneracies between CO2 concentration and surface pressure further complicated the 
interpretation of the calculated spectra, hence the determination of atmospheric conditions. Still, inferring approximative CO2 
concentrations and surface pressures would be possible. 

In practice, detecting atmospheric signals is challenging since calculated SNR values are well below unity in most of the cases. 
^ ' The SNR for a single transit was only barely larger than unity in some near-IR bands for transmission spectroscopy. 

Most interestingly, the false-positive detection of biomarker candidates such as methane and ozone could be possible in low res- 
r***" . olution spectra due to the presence of CO2 absorption bands which overlap with biomarker spectral bands. This can be avoided 

' however by observing all main CO2 IR bands instead of concentrating on, e.g., the 4.3 or 15/jm bands only. Furthermore, a 

C*") ' masking of ozone signatures by C0 2 absorption bands is shown to be possible. Simulations imply that such a false-negative 

qq , detection of ozone would be possible even for rather large ozone concentrations of up to 10~ 5 . 

Key words. Planets and satellites: atmospheres, Stars: planetary systems, Stars: individual: Gliese 581, Planets and satellites: 
individual: Gliese 581 d 
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1. Introduction dWinn et al.ll201 lllDemorv et al.ll201 lh 

Currently, more than 500 extrasolar planets and planet can- The unique geometrical orientation of transiting planets of- 

didates are known. Over 100 of these planets transit their fers the opportunity for the spectral characterization of their 

central star. Ten of the transiting planets are so-called atmospheres. Transmission spectroscopy during the primary 

Sup er-Earths with ma sses below 10 Earth masses: C0R0T - transit favors ultraviolet (UV), visible and near-infrared (IR) 

7 b dLeger et al.ll2009h . GJ 1214b dCharbonneau et al.ll2009l) . wavelengths, since the stellar signal is stronger towards 
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shorter wavelengths. Emission spectroscopy during the sec- 
ondary eclipse is easier in the mid-IR since the planet-star 
flux ratio is higher in this wavelength regime. Thus, both 
methods are complementary. 

Emission and transmission observations have been 
performed for almost 30 extrasolar planets to date. 
Thermal emission of radiation has been detected (e.g., 
Deming et al.l 120051 [Harrington et ail 120061 iKnutson et al] 
20071 lAlonso et al.ll2009l) . Furthermore, the chemical com- 
position of the atmospheres and exospheres of some of 
these planets has been determined. Atoms (H, C, O, Na, 
Fe, Mg) and molecules (CO, CQ 2 , CH 4 , H 2 Q) have been 
detec t ed (e. g.. | Charbonneau et al.l 120021 IVidal-Madiar et al. 



2004[ iTinettiet al.l 12007 
IT 



2009 



Stevenson et al. 



Grillmair et al. 2008, Swain et al 
20101 IMadhusudhan et al.l 120111) . 



iMadhusudhan et al.l 
Although there is some discussion on- going as to whether 
some of these observations are valid dGibson et al. l l20TTl 
iMandell et al.l 1201 ll) . atmospheric characterization of ex- 
oplanets is indeed feasible with current instrumentation. 
For two transiting Super-Earths, spectroscop ic observations 
have been already pe rformed (CoRoT-7 b, Guenther et al.l 
201 ll and GJ 12 1 4 b. iBean et al.1 l20ld iDesert et al.l 1201 ll 
Croll et al l I20TTI ICrossfield et all l201lb . For CoRoT-7 b, 



upper limits on the extension of the exosphere have been 
obtained. The observations of GJ 1214b narrowed down the 
range of possible atmospheric scenarios, favoring either a 
cloud-free water vapor atmosphere or a cloudy hydrogen- 
dominated, methane-depleted atmosphere. These two planets 
are however far too hot to be considered habitable in the 
classical sense of life as we know it on Earth. 
Nevertheless, studies of potential atmospheric signatures 
of terrestrial habitable planets have been performed, in 
order to predict si gnal strengths and assess observa- 
tion s trategies (e.g., Pes Marais et al. 20021 Segura et al.l 
onnil tehr-onr-air-v, at oi I 1^0061 " S ' ^™^Twnu 



2003, Ehrenreich et al 
Miller-Ricci et al.l 12009 
201 lllRauer et alj|201 lb 



Kaltenegger & Traub 2009, 
Deming et al.l l2009l iBelu et afl 



The extrasolar planet GL581d dUdrv et alJ 120071 
iMavor etaL 20091) is the fi rst pot en tially habitable Super- 
Earth (IWordsworth etall 120 1 Obi Ivon Paris et al.l 1201 



Hu & Dindl201 ll IKaltenegger et aljfeoi ll IWordsworth et al 



20101 
i et al.1 



20ill) . The orbital inclination of GL581d has been shown 
to lie in a range between 40-88°, based on photom etric 



constraints for GL581b dLopez-Morales et al] 120061) and 
dynam ical simulations of the whole system dMavor et al] 
120091) . Hence, GL581d is unlikely to transit its central 
star. However, we used GL581 d as an analogue of similar, 
transiting systems which are anticipated to be found in the 
near future. In this study, we illustrate the possible spec- 
troscopic characterization of potentially habitable planets 
based on a wide r ange of atmos p heric model scenarios 
of GL581d from Ivon Paris et all (1201 Oh . following the 
strategy outlined below. Synthetic spectra of some specific, 
potentially habitable GL581d model atmospheres have 



already been presented by IKaltenegger et al] d201 ll) and 
IWordsworth et al] d201 ll) . However, they did not discuss the 
potential detectability (i.e., signal-to-noise ratios) of spectral 
features or the potentially possible detailed characterization 
of their model atmospheres. 

The paper is organized as follows: Section|2]briefly outlines 
the mentioned observing strategy. Section [3] presents atmo- 
spheric scenarios and models used. Results will be shown in 
Sect. [4] Conclusions are given in Sect. [5] 

2. Potential atmospheric characterization of 
terrestrial exoplanets 

2.1. Observations 

The observable quantity during transit measurements is the 
wavelength-dependent planet-to-star contrast ratio, i.e. the 
transit or eclipse depth. Based on the knowledge of the stellar 
radius, transmission spectra thus measure the apparent plan- 
etary radius as a function of wavelength. However, this de- 
pends critically on the accurate characterization of the cen- 
tral star, which is the main source of uncertainty for de- 
rived planetary pr operties (see, e.g., the case of GJ 1214b, 
ICarter et al.ll201l1) . From the stellar properties such as spec- 
trum and radius, and adopting a baseline value for the plan- 
etary radius, emission spectra during secondary eclipse can 
be translated into brightness temperature spectra. The bright- 
ness temperature spectra are particularly illustrative since the 
brightness temperature is the apparent atmospheric temper- 
ature at a given wavelength. For optically thin spectral win- 
dows, the brightness temperature corresponds to the surface 
temperature. 

Note that in the following, we assume that the stellar proper- 
ties as well as the geometric planet radius are known exactly. 

2.2. Atmospheric characterization 

An observation strategy would aim at (1) establishing the ex- 
istence of an atmosphere, (2) determining the major atmo- 
spheric constituents along with radiative trace gases, (3) char- 
acterizing surface conditions, hence assessing habitability 
and ultimately, (4) searching for atmospheric species which 
would indicate the presence of life on a planet. Usually, 
however, the existence of an atmosphere and its composi- 
tion can only be established via atmospheric and spectral 
mod eling with subsequent compari son to the obtained data 



(e.g..lMadhusudhan & Seagel2009llMiller-Ricci et ai]l2009l 



iMiller-Ricci & Fortnevll2010l) 
An additional challenge is the possible false-positive or 
false-negative identification of so-called b iomarkers (e . g., O3 , 
C H4, N?Q) which was discuss ed, e.g., bv lSelsis etal1(l2002l) 
or lSchindler & Kastind (|2000). Biomarkers are atmospheric 
species assumed to be indicative of the presence of a bio- 
sphere on the planet. On Earth, N2O, for example, is be- 
lieved to be almost exclusively produced by denitrifying bac- 
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teria, and O3 is a photochemical product of oxygen which 
itself originates mainly from photosynthetic organisms. Such 
biomarkers are detectable in the Earth's atmosphere due to 
absorption bands (ozone: 9.6 /mi, nitrous oxide: 7.8 and 4.5 
jum, methane: 7.7 and 3.3 fim). 

It has to be distinguished between the false-positive or false- 
negative detection of biomarker species and the false-positive 
interpretation of detected biomarkers as a sign for life. This 
is due t o the fact that abiotic formation of ozone is possi- 
ble (e .g. JSegura et al.ll2007llDomagal-Goldman & Meadows! 
12010) even though the magnitude of the effect is de- 
bated dSelsis et alj|2002l) . Note, for exampl e, that ozone has 
been found in the Martian atmosphere (e.g. lYung & deMorel 
[1991 . Therefore, the detection of ozone alone could still be 
a false-positive detection of life. 

In terms of detecting biomarkers and life, the triple signa- 
ture O3, CO2 and H2O is a possibility to avoid f alse-positive 
detect i ons of biospheres, as proposed by, e.g., ISelsis et ail 
d2002h . lSagan et all (Il993h proposed 2 (or its tracer O3) and 
CH4 as combined biomarkers. 



Model CO2 concentrations were chosen to be consistent 
with CO2 concentrations on present Venus and Mars (95%), 
present Earth (3.55TC T 4 ) as well as assumed scenarios of the 
early Earth (5%, e.g., Kasting i987h . Surface pressures were 
chosen such that model column densities were comparable to 
scenarios of early Earth or early Mar s in the literature (e.g., 
iGoldblatt et al. 20091 iTian et aill20 1(H) . Three different atmo- 
spheric scenarios were considered, defined by the CO2 con- 
centration: low (355 ppm volume mixing ratio, vmr), medium 
(5% vmr) and high (95% vmr) CO2 cases. The model atmo- 
spheres further contained water and molecular nitrogen as a 
filling gas, i.e. the model atmospheres are C02-H 2 0-N2 mix- 
tures. In each of these cases, the surface pressure was varied 
from 1 to 20 bar. For the 1 and 20 bar cases, respectively, 
spectra are presented below. The surface albedo was kept 
constant at A SUI f =0. 1 3 for all runs which is the measured sur- 
face albedo of Earth dRossow & Schifferf[l999h . In doing so, 
clouds were explicitly excluded in the climate calculations of 
von Paris et alJ d2010h . Tabled] summarizes the atmospheric 
scenarios considered in this work. 



3. Models and atmospheric scenarios 

3.1. Atmosphere model 

The spectra shown here have been calculat ed using atmo- 
spheri c profiles from scenarios summarized in lvon Paris et alj 
d2010h . Atmospheric profiles (temperature, pressure, water) 
were calculated using a ID radiative-convecti ve model. The 
model is originally based on the model of iKasting et al.l 
d 1984ft . It solves the radiative transfer equation to calcu- 
late temperature profiles in the stratosphere. The stellar flux 
is treated in 38 spectral intervals of varying width rang- 
ing from 0.237 to 4.545 yu m. The radiative t ransfer scheme 
uses a 5-2-stream method dToon etal.lll989h to incorporate 



Rayleigh scattering by N2, CO2 and H2O ( Vardavas & Carverl 
11984 Ivon Paris et alj2010b . Additionally, absorption by H 2 
and CO? is taken into account in the visible and near-IR 
dPavlov et al.ll2000T) . In the troposphere, the atmosphere is as- 
sumed to be convective. Hence, tempera t ure pr ofile s are as- 
sumed to be adiabatic, based on iKastingl d 1 9881) and lKastingI 
d!991l) . The water p rofile is calculated using a fixed rela- 
tive humidity profile dManabe & Wetheraldll 19671) Mor e de- 
tails on the model ar e given in Ivon Paris et al. r d2008l) and 
Ivon Paris et"ai1d2010l) . 



3.2. Atmospheric scenarios 

The model scenarios used the orbital d istance of 0.22 AU 
and eccentricity of e=0.38 for GL 58 1 d dMavor et al.ll2009h . 
The stella r input spectrum was ba sed on a synthetic NextGen 
spectrum dHauschildt et al.lll999l) in the visible and near-IR 
merged with measurements by the Internati onal Ultraviolet 
Explo rer satellite in the U V, as described by Ivon Paris et al] 
(f2010h . 



Table 1. Atmospheric scenarios for GL 58 1 d considered here 



Scenario 



P [bar] 



I0W-CO2 

medium-C02 

high-CQ 2 



1,20 
1,20 
1,20 



CO2 vmr 
3.55 • 10" 4 

0.05 

0.95 



The model calculations of von P aris et a I] (I2010h suggested 
that four scenarios with high CO2 partial pressures (high- 
C0 2 5, 10 and 20 bar as well as medium-C02 20 bar) 
were habitable, with surface temperatures above 273 K, 
i.e. above the freezing point of water. These results were 
in broad agreement with res ults from | Wordsworth et al.l 
d2010bl) . lHu & Dingl (1201 ll) and lKaltenegger et all (1201 ll) . 



3.3. Computation of spectra and signal-to-noise ratios 

The line-by-line code MIRART-SQuIRRL 
dSchreier & Bottgerl 120031) was used for the calculation 
of high-resolution synthetic planetary emission, bright- 
ness temperature and transmission spectra. Note that the 
collision-induced absorption (CIA) of CO2 is not included 
here, although the CIA is thought to be important in some 
spectral regions around 7 jum and longwards of 20 fim (e.g., 
IWordsworth et al]|2010al) . 

Furthermore, the output of the stellar radiative transfer 
code of the climate model (see above) was used to produce 
reflection spectra /«: 
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Fig. 1. Spectra of the 20 bar high-C02 case. Left: Secondary eclipse spectrum. Inlet shows the corresponding temperature- 
pressure profile. Right: Transit depth spectrum. 



(1) 



where Is is the stellar spectrum, As the spectral albedo from 
the stellar radiative transfer code, R s the stellar radius and 
d the orbital distance. This approach of us ing the climate 
model output to construct spectra is based on Kitzmann et al 



(2011 



2011 



). A similar method was used by IWordsworth et al 



) to calculate synthetic broadband emission spectra of 
GL 581 d model atmospheres. 

These spectra were then used to calculate contrast spectra 
(i.e., emission + reflection) as well as spectra of effective 
tangent height of the atmosphere . Calculations in our work 
closely follow lRauer et al ] (1201 11) where more details can be 
found. The spectra were calculated on an equidistant spectral 
grid (in wavenumber), hence the spectral resolution R varies 
between R ^2-10T0 3 , depending on wavelength. 
For a detection of a spectroscopic feature, the relevant quan- 
tity is the signal-to-noise ratio (SNR). In contrast to the 
SNR calculations by, e.g., IKaltenegger & Traubl d2009l) or 
iRauer et all d201 ll) . we take into account not only the stel- 
lar photon noise, but also the thermal emission of the tele- 
scope, the zodiacal emission and the dark noise. The spec- 
trum of GL 581 used in the SNR calculations is the syn- 
thetic spectrum described above. We base our telescope pa- 
rameters on the James Webb Space Telescope (JWST), as- 
suming a 6.5 m aperture an d a detection efficiency of 0.15 
(IKaltenegger & Traubll2009l) . More details on the noise esti- 
mates can be found in the Appendix [A] The fictitious transit 
duration of GL581 d, hence the assumed integration time, is 
calculated to be 4.15 hours. SNRs are calculated for a spec- 



tral resolution of R-IQ, which is a reasonable value in the 
context of exoplanet characterization. 



4. Results and discussion 

4.1. Presence of an atmosphere 

The left part of Fig. Q] shows the contrast spectrum of the 
high-C02 20 bar case. It is clearly seen that the reflection 
component dominates the near-IR up to about 2.5 fim and in 
the window between 3-4 jum . For longer wavelengths, the 
emission of the planet is the main component. The broad wa- 
ter and CO2 absorption bands are clearly seen in the spec- 
trum. Hence, the existence of an atmosphere as well as the 
presence of water and CO2 could be inferred. Interestingly, 
the planet-star contrast is rather low, even though GL 581 
is an M-type star and GL581 d a Super-Earth. The contrast 
reaches about 4-10^ 5 in the mid-IR which is about an order 
of magnitude higher than the contrast between Earth and the 
Sun. However, it is about 100 times lower than corresponding 
values for hot Jupiters. The right part of Fig.[T]shows the syn- 
thetic transit depth spectrum for the 20 bar high-C02 case. It 
can be seen that due to the presence of large amounts of water 
and CO2, the planet appears larger than its geometric radius 
at all IR wavelengths. Hence, the presence of an atmosphere 
could also be clearly inferred from transmission spectra due 
to the wavelength-dependent apparent radius. 
Figure [2] also shows contrast and transit spectra, but for the 
1 bar I0W-CO2 case. The spectra are relatively flat, except 
for the strong CO2 fundamental bands. In these bands, the 
presence of an atmosphere could be inferred, as for the 20 
bar high-C02 case. Since the atmosphere is very dry (partial 
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Fig. 2. Same as Fig.Q] but for the 1 bar I0W-CO2 case. 




Fig. 3. Brightness temperature spectra (based on planetary emission only) of the 1 and 20 bar scenarios. Left: L0W-CO2 case. 
Right: High-C02 case. Surface temperatures are indicated by horizontal dashed lines. 



pressure of water less than 10~ 5 bar), the water bands are dif- 
ficult to discern in the emission spectrum. In contrast to the 
high-C02 20 bar case, the secondary eclipse spectrum of the 
I0W-CO2 1 bar case is dominated by the reflection component 
up to 4.5 jum, i.e. the limit of the stellar radiative transfer code 
of the atmospheric model. 



4.2. Atmospheric characterization 

After securely detecting the atmosphere, the next step would 
be its characterization (i.e., composition, surface pressure) 
and, from there, assessing the surface conditions, hence po- 
tential habitability. 

4.2.1. Surface temperature 

Figure [3] shows the brightness temperature spectra of the 1 
and 20 bar scenarios of the I0W-CO2 (left) and the high- 
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Wavelength [/im] 




10 15 
Wavelength [firn] 



Fig. 4. Comparison of secondary eclipse spectra of different scenarios: Pressure effect. Scenarios as indicated. Inlet shows the 
corresponding temperature-pressure profiles. 




10 15 
Wavelength [fJ.m] 




Fig. 5. Transmission spectra: Pressure effect. L0W-CO2 (left) and high-CC>2 (right). 



CO2 case (right). Note that brightness temperature spectra are 
based only on the emission spectrum since brightness tem- 
peratures calculated from the contrast due to the reflection 
spectrum would yield values up to about 700 K (near 1 fim). 
Therefore, characterization of atmospheric temperatures of 
terrestrial planets is not possible in the near-IR up to about 
4-5 /im. 

The left part of Fig. [3] shows that for the I0W-CO2 scenarios, 
the surface temperature could be inferred from the brightness 
temperature spectra since the atmosphere is transparent ex- 



cept in the CO2 fundamental bands. By contrast, as can be 
seen in the right part of Fig. [3] the difference between the 
brightness temperature and the surface temperature is always 
non-zero in the high-CC<2 20 bar case. This means that the 
emission spectrum does not allow for a determination of the 
surface temperature, hence to assess potential habitability. 
The reason for this is that the atmosphere is optically thick 
for thermal radiation du e to the large amount s of CO2 and 
water in the atmosphere dvon Paris et alJl2010h . In the 1 bar 
high-CC<2 case, some spectral windows would still allow for 
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the determination of the surface temperature which is mostly 
due to the fact that the atmosphere is much drier than in the 
high-C02 20 bar case. 

4.2.2. Surface pressure 

In Fig. [4] the 1 and 20 bar runs with high and low CO2 con- 
centrations are compared to each other to illustrate how sur- 
face pressures could be inferred from secondary eclipse spec- 
tra. In the I0W-CO2 case (left), the main effect can be seen in 
the 15 /urn band, which is considerably broader for the 20 bar 
run than for the 1 bar run. This is simply due to the fact that 
the line center becomes optically thick at pressures of about 
100 mbar, whereas the line wings are transparent up to pres- 
sures of the order of 5-10 bar. Also, the 2.7 and 4.3 fim bands 
are much deeper in the 20 bar scenario compared to the 1 
bar scenario. In the high-C02 case (right), the pressure effect 
in the 15 fim CO2 band is much less pronounced. Line wings 
are already saturated at pressures well below 1 bar. Due to the 
large difference in surface temperature (see inlet), some spec- 
tral regions (e.g., windows near 8-9 jim or around 1 1 ^m) dif- 
fer considerably. Furthermore, in the near-IR, the reflection 
spectra show large differences, owing to the strong decrease 
of the spectral albedo with surface pressure. This demon- 
strates the use of near-IR secondary eclipse measurements for 
determining atmospheric characteristics besides temperature. 
Figure|5]shows the effect of changing surface pressure on the 
transmission spectrum for the low and the high-CCh 1 and 
20 bar cases. The spectra show significant differences. For 
example, in the 9.5 fim CO2 band, tangent heights differ by 
about a factor of 2 to 3. However, in terms of absolute height, 
this amounts to 5-10 km at most. Note that the CO2 bands at 
7-10//m are visible already in the spectra of the high-pressure 
I0W-CO2 scenarios, contrary to the emission spectra where 
these bands did not appear (see Fig.|4|. 
These results imply that transmission spectra are in general 
more sensitive to surface pressure than emission spectra for 
the cases studied. This is due to two main reasons. Firstly, the 
tangent height to first order depends on the atmospheric scale 
height, H ~ T/m a (T is the surface temperature and m a the 
mean molecular weight of the atmosphere). Thus, for higher 
surface pressures, and corresponding higher surface temper- 
atures, scale heights are larger. Secondly, for higher surface 
pressures, atmospheres extend further out to space. For exam- 
ple, the 20 bar high-C02 case has its model lid (correspond- 
ing to a pressure of 6.6T0 -5 bar) at 20 km, whereas the 1 bar 
case has its model lid at 13 km altitude, which corresponds 
roughly to 3 scale heights difference. 

4.2.3. Surface albedo 

The reflectivity of the surface, the surface albedo, is an im- 
portant parameter to distinguish between different types of 
surfaces, such as oceans or ice. Reflection spectra could in 
principle reveal information about the value of the surface 



albedo. 

However, for optically thick atmospheres, as already dis- 
cussed above, the reflection spectrum does not yield infor- 
mation on the surface albedo in the high-CCh 20 bar case. 
Without the influence of the atmosphere, the reflection con- 
trast would be 7.1 -10 -9 (with a fixed surface albedo of 
^surf=0.13, see Sect. I3.21) . As can be seen from Fig. Q] the 
reflection contrast is far lower than this value, indicating that 
even cloud-free atmospheres could inhibit surface character- 
ization via reflection spectra. The same is true for the high- 
CO2 1 bar and the I0W-CO2 20 bar cases (not shown). Only in 
the I0W-CO2 1 bar case, the reflection spectrum would allow 
to infer the surface albedo (see Fig. [2J. 




Wavelength [/a 



Fig. 6. Comparison of secondary eclipse spectra of differ- 
ent scenarios: Concentration effect. Inlet shows the used 
temperature-pressure profiles. 



4.2.4. Atmospheric composition 

Fig.[6]shows secondary eclipse spectra at equal surface pres- 
sures (20 bar), but for different CO2 concentrations (high, 
medium, low). It is rather difficult to distinguish between 
the medium and high-C02 scenario, except for some narrow 
spectral windows. When comparing Fig. [4] with Fig. [6] it is 
therefore difficult to decide whether the shape of the spectrum 
is actually due to a pressure difference at high CO2 concentra- 
tion or a concentration difference at high surface pressures. It 
is, however, possible to distinguish between high or low CO2 
concentrations by the presence of many weak CO2 bands in 
the high-C02 case, e.g. at around 7, 9 and 10/jm. These bands 
do not appear in the spectrum unless CO2 concentrations ex- 
ceed several percent. Additionally, the presence of water can 
be inferred from the strong 6.3 /im fundamental band and the 
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signatures of the rotation band in the mid- to far-IR. Thus, 
dry and moist atmospheres could be distinguished. 
In order to show the effect of changing CO2 concentration 
on transmission spectra, Fig.|7]compares the 20 bar scenarios 
for low, medium and high CO2 concentrations. Interestingly, 
the medium-CC>2 case shows a more pronounced 15 /mi band 
of CO2 despite the fact that the high-CC>2 scenario is signifi- 
cantly warmer at the surface (by about ~ 40 K). This is due to 
the lower mean atmospheric molecular mass (hence, higher 
scale height) which is ~ 29 g moF 1 in the medium-CCb case 
and ~ 43 g moF 1 in the high-CC>2 case. The weak bands 
around 7 and 10 /mi do not differ by much. Still, in the wa- 
ter rotation bands (longwards of 20 /mi) and the strong CO2 
fundamentals, rather large differences in the spectrum can be 
seen. This indicates that transmission spectroscopy could be 
able to distinguish between different CO2 concentrations. 




Wavelength [ftm] 



Fig. 7. Transmission spectra: Concentration effect. 

The same is true for the comparison of low and high-CC>2 
scenarios. Here, however, the differences in the spectrum are 
visible in the weak bands of CO2 (e.g., around 10 //m) rather 
than in the strong fundamentals (4.3 and 15 /mi). 

4.2.5. Surface conditions 

In summary, the characterization of surface conditions on 
GL581 d by secondary eclipse spectroscopy would be rather 
difficult, even if the planet were to be transiting. It would 
be possible to approximately constrain CO2 and water con- 
centrations in the atmosphere, however, constraining sur- 
face pressures or surface temperatures, hence assess hab- 
itability, is complicated by degeneracies, as stated above. 
Overall, results imply that it is easier to characterize the 
atmospheric scenarios of GL581d with transmission spec- 



troscopy than with emission spectroscopy, especially with 
respect to the habitable scenarios with massive CO2 at- 
mospheres. However, surface pressures cannot be inferred 
directly with transmission spectra since effective tangent 
heights are always of the order of a few kilometres. Still, 
when combining several spectral bands and using both trans- 
mission and emission spectra, it may be possible to constrain 
the surface pressure as well as CO2 concentrations and the 
presence of water. Thus, in principle, through atmospheric 
modeling surface conditions could be assessed. 

4.3. Biomarkers 

Although the model scenarios used in this work only con- 
sider H2O-CO2-N2 atmospheres, we can compare the com- 
puted spectra with spectral signatures of modern Earth to dis- 
cuss the potential for false-positive or false-negative identifi- 
cations of biomarkers. 

4.3.1 . False-positive detections of biomarkers 

As an illustration, Fig. [8] shows high-resolution and binned 
(with R-10) emission spectra of GL581d in comparison 
to modern Earth. T he modern Earth spectrum (taken from 
iGrenfell et"ai]|2011l) is shown in the upper left, the GL 581 d 
scenarios from this work are (clockwise) high-C02 20 bar, 
high-C02 1 bar and the medium-C02 20 bar case. The spectra 
are centered around 9.6 /mi which is the position of the main 
absorption band of ozone. This absorption feature is close to 
an absorption band of CO2 at around 9.5 /mi which is clearly 
seen in all GL 58 1 d spectra in Fig. [8] 

Additionally indicated is the position of the 7.7 //m band 
of methane (overlapping the 7.8 /mi band of nitrous oxide) 
which again is close to absorption bands of CO2 at 7.5 and 
8 /mi. The methane band is hardly discernible in the Earth 
spectrum, at the considered resolution of R-10. However, 
in the GL581 d scenarios shown in Fig. [8j an absorption is 
seen, which in this case is due to CO2. These CO2 bands at 
biomarker positions (7.5, 8, 9.5 /urn) are also present in the 
transmission spectra (e.g., Fig. |7), implying that transmis- 
sion spectra also suffer from the possibility of false-positive 
biomarker detections. 

At this low spectral resolution, the spectral features around 
9.6 //m, due to CO2 in the GL581 d cases or due to ozone 
in the Earth case, look similar. Thus, the CO2 bands could 
be mistaken for actual biomarker signals, at the expected low 
spectral resolution used for exoplanet characterization. This 
is the case for the emission spectra of medium and high-C02 
scenarios, whereas in the I0W-CO2 case (not shown in this 
Fig.), the respective CO2 bands are too weak to produce false 
detections. In the transmission spectra (e.g., Fig. [7]), these 
bands also appeared for high-pressure, I0W-CO2 scenarios 
which implies that false-positives are possible in transmis- 
sion spectra even for these scenarios. 
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Fig. 8. Possibility of a fals e-positive biomarker detection for different scenarios, as indicated. Modern Earth spectrum taken 
from lGrenfelletail j20TTb . Planck curves of corresponding surface temperatures as dashed lines. Smoothed (red) and binned 
spectra (blue) are for a spectral resolution of 10. 



One possible way of avoiding such false-positive detections 
in the cases presented here (i.e., CCh-rich atmospheres) is 
to exploit the double nature of the CO2 bands around 7 and 
10 //m. If spectral observations are performed, e.g., at 9.5 
and 10.5 /vm, and both filters show a deep absorption, then 
the spectral signatures are most likely due to CO2. Hence, a 
possible spectral characterization with respect to biomarkers 
should be done in all main IR CO2 bands (2-15 /mi) in order 
to avoid false-positive detections. 



4.3.2. False-negative detection of biomarkers 

In order to investigate the possibility of false-negative detec- 
tions of ozone, we calculated model spectra of the high-C02 



20 bar case and of the medium-CCh 20 bar case where we 
additionally introduced artificial profiles of ozone in the line- 
by-line spectral calculations. These ozone profiles have not 
been implemented in the ID climate calculations. However, 
they are not expected to strongl y influence the tempera ture 
structure of the atmosphere (e.g. jKaltenegger et al.l201 it) be- 
cause of the lack of UV radiation emitted by GL 581 in the 
wavelength regime where ozone heating occurs. 
In a first attempt, we used isoprofiles of ozone in the calcula- 
tion of the spectra. Chosen concentrations were 10~ 8 , 10 -7 , 
10~ 6 and 10~ 5 , respectively. This choice of concentrations 
broadly covers the range of ozone concentrations found in 
the present Earth atmosphere (10~ 8 near the surface, about to 
almost 10~ 5 in the mid-stratosphere). 

Results of this sensitivity analysis are shown in Fig.|9]for the 
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Fig. 9. Effect of O3 isoprofiles on emission spectra. Constant O3 concentrations as indicated. Left: High-CC<2 20 bar, right: 
Medium-CCh 20 bar. Horizontal lines represent binned spectra for a spectral resolution of 10. 



high-C02 and medium-C02 20 bar cases. Clearly, a detection 
of ozone at low concentrations of 10~ 8 or 10~ 7 would not be 
possible at low spectral resolution. For the high-CCh 20 bar 
case, even a concentration of 10~ 6 would be very challeng- 
ing to detect. In both cases shown here, however, it would 
be possible to infer ozone levels of 10~ 5 since the calculated 
intensity drops by about a factor of 3-8 in the center of the 
ozone fundamental band. 

In reality, ozone profiles will most likely not be in the form 
of isoprofiles. The production of ozone proceeds through the 
three -body reaction 

2 + + M->0 3 +M (2) 

where M is, on Earth, typically nitrogen. The production of 
the atomic oxygen needed in eq. [2] requires photolysis of 
molecular oxygen, carbon dioxide or water. Hence, there is a 
trade-off between the necessary UV radiation for photolysis 
and the density to enable the three-body production reaction, 
which is why photochemical models generally predict a dis- 
tinct maximum of atmospheric ozone concentrations around 
the 1-10 millibar pressure range. 

The actual location of this ozone maximum depends on the 
stellar UV radiation fi eld, the ozone and the oxygen con- 
tent of the atmosph e re (ISelsis et al.ll2002llSegura et alj | 2003l 



Segura et al 



Grenfell et al 



20071 iDomagal-Goldman & Meadows! bold 



2011). Therefore, we inserted artificial ozone 



profiles Co 3 as a function of pressure p based on a Gaussian 
profile 



Co, (p) = C max ■ exp 



(*< )' 
0.5 



2\ 



(3) 




0.010 0.100 1.000 

3 concentration [ppm] 



10.000 100.000 



Fig. 10. Mo d ern E arth mean ozone profile from 
iGrenfell et al] (1201 ll) and approximation (in red) with 
eq. fusing C max =9ppm and p„,-6 mbar. 



where C max is the maximum ozone concentration reached at 
pressure p m . Fig. [10] shows the approximation of an Earth 
ozone profile with eq. [3] The agreement near the maximum 
of the ozone layer is rather good. 

FigureQT]shows the effect of varying C max on emission spec- 
tra at constant p m of 10 mbar. Shown are the spectra for the 
high-C02 and medium-CCh 20 bar cases. In Fig. [12] the ef- 
fect of varying p m on the spectrum at constant C max of 10 
ppm is illustrated for the high-C02 20 bar case. 
Note that tropospheric concentration maxima, i.e. at pres- 
sures higher than 100 mbar, are usually not expected from 
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Fig. 11. Effect of O3 profiles on emission spectra. Parameters for eq. [3] C max as indicated, /;>,„= lOmbar. Left: High-C02 20 
bar, right: Medium-CCh 20 bar. Horizontal lines represent binned spectra for a spectral resolution of 10. 



atmospheric chemistry modeling. Rather, maxima are found 
in the lower to upper stratosphere. Hence the scenarios with 
p m of 100 or 1,000 mbar shown in Fig. [12] are most likely 
overestimating the ozone column. 

From Fig. QT| it is evident that it is very difficult to infer the 
presence of ozone even at high concentrations of about 10 
ppm. Only for the profile with C max =100ppm, it would be 
possible to detect an additional absorption besides the CO2 
band. 

Figure [T2l shows that for a value of /? m = 100mbar or higher, it 




would be possible to detect ozone in the spectrum. However, 8 8 5 9 9 5 10 10 5 

r r Wavelength [fim] 

as stated above, such a tropospheric maximum of ozone is 
not very probable. 



Fig. 12. Effect of maximum of O3 profiles on emission spec- 
tra for the high-C02 20 bar case. Parameters for eq. [3] 
Cmax=10ppm, p m in mbar as indicated. Horizontal lines rep- 
resent binned spectra for a spectral resolution of 10. 



These sensitivity studies illustrate that an inferred absence of 
ozone could be possible for the atmospheric scenarios con- 
sidered in this work, even if ozone would be present in large 
concentrations of the order of 10 6 -10 5 . Hence, also the pos- 
sibility of false-negatives of ozone needs to be taken into ac- 
count when investigating COa-rich, potentially habitable at- 
mospheric scenarios. 



4.4. Detectability 



The main spectral bands investigated for detectability are the 
2.0, 2.7, 4.3, 7.7, 9.5 and 15 fim CO2 absorption bands as well 
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Table 2. SNR values for secondary eclipse spectroscopy. Telescope configuration JWST, integration time 4. 15 hours, spectral 
resolution R=10. 



Scenario 


8.5 yum 


9.5 yum 


11.5 yum 


15 ,um 


20 yum 


I0W-CO2 1 bar 


0.050 


0.091 


0.19 


0.10 


0.23 


I0W-CO2 20 bar 


0.062 


0.10 


0.22 


0.088 


0.25 


medium-C02 1 bar 


0.064 


0.10 


0.22 


0.091 


0.26 


medium-C02 20 bar 


0.56 


0.28 


0.51 


0.089 


0.34 


high-C0 2 1 bar 


0.12 


0.15 


0.37 


0.10 


0.34 


high-CQ 2 20 bar 


0.77 


0.22 


0.74 


0.10 


0.48 



Table 3. SNR values for transmission spectroscopy. Telescope configuration JWST, integration time 4.15 hours, spectral 
resolution R=10. 



Scenario 


2.0 yi/m 


2.7 yum 


4.3 yum 


6.3 yum 


7.7 /im 


9.5 yum 


15 yum 


I0W-CO2 1 bar 


0.41 


0.65 


0.83 


0.091 


0.027 


8.4 10- a 


0.21 


I0W-CO2 20 bar 


1.1 


1.2 


1.4 


0.24 


0.22 


0.12 


0.33 


medium-CCb 1 bar 


1.0 


1.1 


1.2 


0.12 


0.23 


0.12 


0.29 


medium-C02 20 bar 


2.1 


1.9 


2.0 


0.46 


0.65 


0.46 


0.44 


high-C0 2 1 bar 


0.96 


0.89 


0.94 


0.10 


0.26 


0.15 


0.22 


high-CQ 2 20 bar 


1.9 


1.6 


1.6 


0.50 


0.65 


0.46 


0.35 



as the 6.3 /urn H2O absorption band in transmission. For sec- 
ondary eclipse spectroscopy, only the 9.5 and 15 /im CO2 and 
the 20 yum H2O bands were considered. Additionally two fil- 
ters outside of broad absorption bands are calculated, namely 
at 8.5 and at 11.5 yum. These two filters offer the possibility 
to characterize surface temperatures in optically thin atmo- 
spheres. 

Furthermore, six representative scenarios are considered, 
namely the 1 and 20 bar runs of the low, medium and high- 
CO2 cases. 

Table [2] summarizes the SNR values for secondary eclipse 
spectroscopy. It is clearl y seen that, even th ough GL 581 is a 
very close star (6.27 pc. iButler et~ai1l2006l) . obtainable SNR 
values are very small, only reaching values of 0.3-0.77. The 
reason for the low SNR values in the mid- to far IR is the 
effect of the zodiacal background chosen here (see Appendix 
lAl which reduces SNR values by about a factor of 5 at 20 yum 
and a factor of 2 at 15 yum compared to the pho ton-limited 
case. A similar result was already found by, e.g., iBelu et all 
d20Tl . 

Table [3] summarizes the SNR values for transmission spec- 
troscopy. SNR values are mostly below unity. The exceptions 
are the 2.0, 2.7 and 4.3 yum near-IR fundamentals of CO2. 
In these bands, characterization of the GL581d scenarios, 
as outlined above, could be feasible. A clear distinction be- 
tween, e.g., high and I0W-CO2 cases is still very difficult, 



given that the SNR values are only marginally larger than 
unity. One strategy to obtain larger SNRs is the co-adding of 
transits. This might result in a significant increase of the SNR 
since the SNR scales with the square root of the number of 
transits. Assuming a 5-year mission lifetime of JWST (3-4 
transits observable per year) would lead to SNRs of about a 
factor of 4 higher compared to the values shown in Tables [2] 
and[3] 

5. Conclusions 

We presented synthetic emission and transmission spectra 
for a wide range of atmospheric scenarios of the potentially 
habitable Super-Earth GL581 d. These spectra were used as 
an example of possible future spectroscopic investigations of 
candidate habitable planets. 

Water and carbon dioxide could be clearly seen in the calcu- 
lated spectra due to prominent absorption bands, indicating 
the presence of an atmosphere. 

The determination of surface temperatures was possible for 
model atmospheres with either low surface pressures or low 
CO2 content. The potentially habitable, C02-rich scenarios 
did not allow for the characterization of surface tempera- 
tures. Thus, their potential habitability could not be assessed 
directly fr om the spectra, in agreem ent with a model spec- 
trum from lKaltenegger et al.l d201 lb . The further determina- 
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tion of atmospheric conditions was complicated by degenera- 
cies between the surface pressure and the CO2 concentration. 
However, when combining observations in several spectral 
bands and using both transmission and emission spectra, in- 
ferring approximative CO2 concentrations and surface pres- 
sures would be possible. 

With currently planned telescope designs such as JWST, SNR 
values for emission and transmission spectroscopy are, how- 
ever, rather low, implying that the detection of an atmo- 
sphere of a "GL581 d-like" transiting planet is challenging. 
Reasonable, single-transit SNR could only be calculated for 
three near-IR CO2 fundamentals at 2.0, 2.7 and 4.3 fim (trans- 
mission spectroscopy). In the mid- to far-IR, thermal and zo- 
diacal background noises inhibit SNRs above unity. 
Results indicate that the search for biomarkers in CCh-rich 
or high-pressure atmospheres would suffer from the possibil- 
ity of false- positive detections, in agre e ment with previou s 
studies (e.g JSchindler & Kastind2000l or lSelsis et al.ll2002l) . 
This is due to absorption bands of CO2 which occur close to 
main biomarker absorption bands. However, if the main CO2 
IR bands were to be observed simultaneously, such false- 
positive detections could possibly be avoided. This will be 
subject of future modeling including Earth-like atmospheric 
composition. Results also imply that CO2 absorption bands 
could mask the spectroscopic features of ozone, hence pro- 
duce false-negative detections. This was shown to be possible 
even if ozone would be present in rather large concentrations 
of up to 1(T 5 . 
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Appendix A: Noise contributions 

The total SNR S p of a planetary feature is calculated with 



S p = 



V total 



(A.l) 



where F p is the planetary signal measured during secondary 
eclipse (emitted and reflected spectrum) or primary eclipse 
(additional transit depth due to atmospheric absorption). cr tota i 
represents the total noise of the measurement. In addition to 
the stellar noise cr s , main contributions to the noise budget 
come from the zodiacal emission of the solar system (denoted 
Czodi), the thermal emission of the mirror and, in the case of 
the JWST, the sun shade (denoted cr therma i). Both components 
need to be taken into account in the SNR calculations since 
their effect on detectability could be sig nificant, especially in 
the mid-IR, lon gwards of 10 //m (e.g.. iDeming et aDl2009l 
lBeluetalJl201lh . In addition, the dark noise <x (/ is considered 
here. Hence, we have 



ctotai _ y2 ■ ( 
The thermal noise was obtained from 



22 2 

5 ^"thermal ^"zodi 



thermal 



q e • £ ■ Bj(T S wst) ■ t, ■ n px ■ S a ■ A e 



A 
R 



(A.2) 



(A.3) 



where e is the emissivity, B,i(7twst) the blackbody emis- 
sion of the telescope and sun shade, t, the integration time 
(here assumed to be the transit duration, i.e. 4.15 hours), n px 
the number of pixels used for the integration on the detec- 
tor, A the considered wavelength and R the spectral resolu- 
tion, Ae the emittin g area, g e =0.15 the quantum efficiency of 
the te lescope (e.g.. iKaltenegger & Traubll2009l iRauer et al.l 
1201 11) and S a the solid angle per pixel. 
We assumed a telescope and sun sha de temperature Tj wst of 
45 K with an emissivity of e=0.15 dBelu et al.ll2011l) and a 
total emitting surface of sun shade and mirror combined of 



A £ =240 m 2 (INella et al.ll2002l) . The number of pixels n px is 



calculated with 



spatial ' "spectral 



R 



(A.4) 



where « spa tiai and n spe ctrai the number of pixels in spatial 
and spectral (per ^m) direction, respectively. We assumed 
"spectral =3 2. 3 3 pixel fxm~ , mi micking the MIR I instrument 
in the 5-11 jum spectral range dBelu et al.ll20TTh . The number 
of pixels in the spatial direction is obtained from 



^spatial 



2.4- 



,1 



df p s 



(A.5) 



with d, the telescope diameter (i.e., 6.5 m in the case of 
JWST) and p s —0.\V per pixel the instrumental pixel scale. 
The number calculated from eq. IA.5l is then rounded up to the 
closest multiple of 4. The solid angle per pixel is calculated 
with (1 sq. deg « sr) 



S A = 



3600 2 • 3282 



(A.6) 



The zodiacal emission Iz (given in W irT 2 pnT 1 sr _1 ) used 
in the noise calculations is taken from o ne example m e asure- 
ment in the ecliptic plane presented by iKelsall et al.l d 1998b 
(their Fig. 9). Note that this choice is a rather pessimistic as- 
sumption since some potential targets will probably be lo- 
cated towards higher ecliptic latitudes, hence the correspond- 
ing zodiacal noise would be somewhat reduced. We calculate 
the zodiacal noise via 



°"zodi - 7 z ■ 'V ' S a ■ Atei t t q e - (A.7) 

with A t ei the telescope area, assuming a circular aperture of 
6.5 m. 
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The dark noise contribution to the noise budget is calculated 
from 

cr 2 d = d c - n px • t, (A. 8) 

wh ere of c =0.03 e pixel 1 s 1 is the dark c urrent (MIRI, see 
e.g. lDeming et al.ll2009l or lBelu et al.ll201 ll) . 
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